Stabilized Galerkin approximation of convection-diffusion-reaction equations: discrete maximum principle and convergence

We analyze a nonlinear shock-capturing scheme for H1-conforming, piecewise-affine finite element approximations of linear elliptic problems. The meshes are assumed to satisfy two standard conditions: a local quasiuniformity property and the Xu–Zikatanov condition ensuring that the stiffness matrix associated with the Poisson equation is an M -matrix. A discrete maximum principle is rigorously established in any space dimension for convectiondiffusion-reaction problems. We prove that the shock-capturing finite element solution converges to that without shock-capturing if the cell Péclet numbers are sufficiently small. Moreover, in the diffusion-dominated regime, the difference between the two finite element solutions super-converges with respect to the actual approximation error. Numerical experiments on test problems with stiff layers confirm the sharpness of the a priori error estimates.